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Even-Odd Oscillation in Conductance of Single-Row Sodium Nanowire 
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We present a first-principles calculation of the electronic conduction properties of single-row 
sodium nanowires suspended between semi-infinite electrodes. The conductance of the nanowire 
is ~ 1 Go (= 2e^//i) and oscillates with a two-atom period as the number of the atoms within the 
nanowire varies. Moreover, we observed bunches of high electron density with a two atom-lengths in 
the channel density distribution. The relation between the period of the conductance oscillation and 
the length of bunches are examined by using simplified models and is found to be largely affected 
by the characteristics of the infinite wire. 

PACS numbers: 73.63.Nm, 68.65.La, 73.21. Hb, 73.23.Ad 



The properties of condensed systems with nanometer 
dimensions have attracted considerable interest in the 
past decade because of their great importance in funda- 
mental physics as well as in nanometer-scale fabrication 
technologies. When the systems' dimensions reach the 
order of the electron mean free path, the electron trans- 
port becomes ballistic and its conductance is quantized 
in the unit of Gq (= 2e'^/h) contrary to the diffusive 
transport in macroscopic systems, where e is the elec- 
tron charge and h is Planck's constant. Since clearly 
showing such unique quantum phenomena, the metallic 
nanowire contacts have been examined in many experi- 
mental and theoretical studies. In experiments, metal- 
lic nanowires are generated using scanning tunneling mi- 
croscopy or mechanically controllable break junctionsjjy. 
Many experiments observed that the nanowires made of 
monovalent atoms such as NaflQ, Cu|l|l,|||, Ag|10 
or Au0, IE S 111 manifest conductance quantization. 
Moreover, Smit et al.^ found in their recent experiment 
that the conductances of Au, Pt and Ir nanowires os- 
cillate with a two-atom period as the number of atoms 
constituting the nanowire varies, which is the so-called 
even-odd oscillation. On the theoretical side, there have 
been several first-principles studies on the electron con- 
duction properties of single-row nanowires, and some 
of them reported that the conductances of the single- 
row nanowires suspended between electrodes exhibit the 
even-odd oscillation[l3 lll|ll|ll|ll|ll|ll|ll|ll|ll 
l20i |. Although these results are consistent with those of 
previous experiments and give us a certain knowledge, 
there remains much to be learned about the origin of the 
conductance oscillation. 

In this Rapid Communication, we carry out a first- 
principles study on the electron conduction properties of 
single-row sodium nanowires for a detailed interpretation 
of the conductance oscillation. We observed that the 
nanowires manifest the even-odd oscillation of conduc- 
tance and the channel charge-density distribution forms 
bunches of high electron density along the nanowire. 
Moreover, the length of the bunches is equal to the pe- 
riod of the oscillation, 2d, where d {d = 7.0 a.u.) is 
the nearest-neighbor interatomic distance in the sodium 
bulk. The relation between the conductance oscillation 



and the emergence of the bunches is addressed by using 
simplified models. 

Our first-principles calculation code is based on the 
real-space finite-difference method [2ll 12^. l23| | . Since this 
approach does not need any basis-function sets, unlike 
the plane-wave approach or linear-combination-atomic- 
orbital (LCAO) method, one can eliminate drawbacks 
due to basis sets, e.g., a strict treatment of the systems 
with the combination of periodic and nonperiodic bound- 
aries is computationally difficult in the plane-wave ap- 
proach and the completeness of the basis set is always a 
concern in the LCAO method. These advantages allow 
us to calculate the conduction properties of nanostruc- 
tures suspended between semi-infinite electrodes with a 
high degree of accuracy. 

We first examine the conduction properties of the 
single-row sodium nanowires sandwiched between semi- 
infinite Na(OOl) crystalline electrodes. Figure ^ shows 
the computational model in which the case of Natom = 3 
is depicted as an example, where Natom is the number 
of atoms in the nanowire. The scattering region consists 
of the single-row sodium wire, the square bases made of 
sodium atoms, and four sodium (001) plane of electrodes. 
The nanowire is connected to the bases on both ends 
and all of them are suspended between the electrodes. 
The distance between the electrode surface and the ba- 
sis, as well as that between the basis and the edge atom 
of the nanowire, is d/^/3 and the distance between adja- 
cent atoms in the nanowire is d. We treat a local pseu- 
dopotential of a sodium 3s electron for the nucleus con- 
structed using the algorithm by Bachelet, Hamann, and 
Schliiter^J. Exchange-correlation effects are treated by 
the local density approximation l25l w ithin the framework 
of the density functional theorv|2q. The central finite- 
difference formula is employed for the derivative aris- 
ing from the kinetic- energy operator in the Kohn-Sham 
eauation[27j. The grid spacing is taken to be 0.81 a.u., 
which corresponds to a cutoff energy of 15 Ry in the 
plane-wave method. The Kohn-Sham effective potential 
is determined self-consistently using the conventional pe- 
riodic supercell which is represented by the rectangle in 
Fig. n In the calculation of the global wave functions for 
infinitely extended states continuing from one electrode 
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side to the other, the overbridging boundary-matching 
method[23 is adopted with the supercell imposed peri- 
odic boundary conditions in the x and y directions and 
the nonperiodic boundary condition in the z direction. 
The electron transmission of the nanowire at the zero 
bias hmit is calculated by the Landauer formula|29j . 
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FIG. 1: Schematic view of single-row sodium nanowire 
[Natom ~ 3) suspended between semi-infinite Na(OOl) crys- 
talline electrodes. The rectangle represents the supercell em- 
ployed to determine the Kohn-Sham effective potential and 
the scattering region to calculate the conduction property. 



Figures El and 121 depict the conductances and the chan- 
nel charge-density distributions of incident electrons with 
the Fermi energy from the left-hand electrode, respec- 
tively. The planes shown in Fig. O are perpendicular to 
the [100] direction and include the nanowire axis. The 
conductances of nanowires are ~ 1 Gg, which are in 
good agreement with other theoretical and experimen- 
tal resultsU |l|li|ll|lllil|ll|ll|ll|ll|ll|23. 
In Fig. 13 one can evidently recognize that the conduc- 
tances of the odd-number nanowires {Natom — 3 and 
5) are higher than those of the even-number nanowires 
{Natom = 2 and 4); the conductance oscillation with a pe- 
riod of two atom-lengths is consistent with that reported 
in the other theoretical studieSjlOt. .IL .12^ .14 .15^ .16i . il7i . 
El in Hi and experiment 0. Additionally, in the case 
of the odd-number nanowires, one observes bunches of 
high electron density with a length of 2d [Figs. 3(b) and 
Old)] which corresponds to the period of the even-odd 
oscillation, and the density distributions of the incident 
electrons are almost symmetric with respect to the (001) 
plane including the center atom of the nanowire. 

Next, we will explore the relationship between the 
even-odd oscillation and the bunches of high electron den- 
sity. For simplicity, we first consider the one-dimensional 
problem of the penetration of a square potential barrier 
as shown in Fig. 4(a), where the nanowire is imitated 
by the potential barrier. An incident electron with an 
energy E has a wave vector of fc = \/2{E — Vq) in the 
barrier region and its transmission T is determined by 
wave-function matching[3l| 
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FIG. 2: Conductances of single-row sodium nanowires. 



where Vq is the height of t he p otential barrier, L is the 
barrier length and k — \J2E. When L is an integral 
multiple of 7r/fc, the transmission shows unity. The elec- 
tron density forms bunches with a length of 7r/fc so that 
it becomes almost symmetric with respect to the center 
of the nanowire. On the other hand, when L is half- 
integral multiple of 7r/fc, the transmission exhibits min- 
imum and the density distribution becomes asymmet- 
ric due to the fractured bunches. More realistically, we 
then consider a three-dimensional prismatic jellium wire 
with side lengths of W connected to semi-infinite elec- 
trodes [Fig. 4(b)]. The potentials for electrons are de- 
fined as zero inside the nanowire and electrodes, whereas 
that in the vacuum regions is infinitely high. It is well 
known that the lowest band energy of the infinite pris- 
matic jellium wire with W in the x and y directions 
\s E = \k1 + {tt/WY, where is the z component 
of the wave vector. When the wire is connected to 
the electrodes, the incident electrons with the energy 
E incoming from the electrode have a wave vector of 
kz — \/2\E — {tt/WY] in the nanowire region and the 
transmission T is determined by 
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kl{p + (3*f 



kl{f3 + (3*Y + \kl-(3-^\^ sm'[kz{L - AL)] • 
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where (3 — — i^^*"(0) with ^™{z) being the channel 
wave function of incident electrons from the left-hand 



(k,+l3){k,-l3') 



Note that 



electrode and AL = tt-^ log 

^''"(2;) includes left-decreasing evanescent waves as well 
as right-propagating waves. Figures \S\ and show the 
transmission as a function of the nanowire length and 
channel charge-density distributions of the incident elec- 
trons from the left-hand electrode, respectively. The pe- 
riod of the oscillation and the length of the bunches do 
not depend on the energy of incident electrons, when they 
are measured in units of Tr/fc^. Furthermore, when L is 
an integral multiple of Tr/fc^, the conductance is 1 Go 
and the channel density distribution is symmetric with 
respect to the center of the nanowire. On the other hand, 
when i is a half- integral multiple of n/kz (Fig. O, the 
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FIG. 3: (color online). Channel charge- density distributions of incident electrons with the Fermi energy from left-hand electrode 
for (a) Natom = 2, (b) Natom ~ 3, (c) Natom ~ 4 and (d) Natom ~ 5. Each contour represents twice or half the density of the 
adjacent contour lines. The lowest contour represents 4.84 x 10"'' electron/bohr^/eV. 



conductance exhibits minimum and the distribution be- 
comes asymmetric. 
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FIG. 5: Conductance of prismatic jellium nanowire as func- 
tion of wire length L. Here, T™i„ is [1 -|- |fc? - /3^p/fc^(/3 -|- 
(a), (b) and (c) are the points where the channel 
charge-density distributions are shown in Fig.|S| 



FIG. 4: Schematic view of (a) potential barrier in one- 
dimensional system and (b) prismatic jellium wire connected 
to semi-infinite jellium electrodes. 



Now, we return to the case of the atomistic sodium 
nanowircs. Since the energy band of the infinite single- 
row sodium wire with the interatomic distance of d 
crosses the Fermi energy at the wave vector kz of 
7r/2ci[33|, the period of the conductance oscillation and 
the length of the bunches become 2d, which justifies the 
above-mentioned explanation. 

Smit et al.^ and Lee et discussed the even- 

odd oscillation behavior using one-dimensional nanowires 
with finite and infinite lengths, and however they did not 
mention the emergence of the bunches of high electron 
density with relation to the conductance oscillation. 



In summary, we have studied the relation between the 
even-odd oscillation of conductance and the bunches of 
high electron density by calculating the conduction prop- 
erties of atomistic sodium nanowires and jellium ones. 
The conductance oscillation can be understood in terms 
of the characteristics of the infinite wire and the conven- 
tional wave-function matching: when the nanowire whose 
length is an integral multiple of Tr/fc^ is connected to the 
electrodes, the conductance is maximum and bunches 
with a length of Tr/fc^ come about so that the channel 
density distribution becomes symmetric. On the other 
hand, in the case where the nanowire length is a half- 
integral multiple of tt/Zcz, the conductance is minimum 
and the distribution becomes asymmetric. The work in 
progress is the investigation of nanowires consisting of 
multivalent atoms. 
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FIG. 6: (color online). Channel charge-density distributions of incident electrons from left-hand electrode at (a), (b) and (c) 
in Fig.|^ Each contour represents twice or half the density of the adjacent contour lines. 
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